FACTORIZED DOMAIN WALL PARTITION FUNCTIONS IN 
TRIGONOMETRIC VERTEX MODELS 



O FODA, M WHEELER AND M ZUPARIC 



Abstract. We obtain factorized domain wall partition functions for two sets 
of trigonometric vertex models: 1. The A'^-state Deguchi-Akutsu models, for 
A^ e {2, 3, 4} (and conjecture the result for all N > 5), and 2. The sZ(r+l|s+l) 
Perk-Schultz models, for {r, s S N}, where (given the symmetries of these 
models) the result is independent of {r, s}. 



0. Introduction 

Domain wall partition functions (DWPF's) were first proposed and eval- 
uated in determinant form, for the spin-i vertex modeQ on a finite square lattice, 
in [1,2]. At the free fcrmion point of the spin-^ model, this determinant is in 
Cauchy form and therefore factorizes. More recently, determinant expressions for 
the DWPF's of spin--^^--!- models and also of level-1 affine so{N) models (for certain 
discrete values of the crossing parameter) were obtained in [3] and in [4], respec- 
tively. 

State variable conjugation. We are interested in models with state variables 
{a}. Each state variable takes discrete integral values, a S {1, • • • , N}. We define 
'state variable conjugation' as replacing each state variable a by {N — ct -I- 1). The 
models mentioned above are invariant under this conjugation. 

The iV-state Deguchi-Akutsu (A^-DA) models, > 2 [5], are models with 
vertex weig hti, that depend on two sets of parameters: 1. Vertical and horizontal 
rapidities, and 2. Vertical and horizontal external field variables. They reduce 
in the limit of no external fields to the spin- ^""'^ models at their respective free 
fermion points. 

Factorized DWPF. In [6], a determinant expression for the DWPF of the 2- 
DA model was obtained using the arguments of [1,2], but only for zero values of the 
rapidities. This determinant is in Cauchy form and therefore factorizes. For general 
values of all parameters, no determinant expression was found, and it was argued on 
general grounds that no such expression exists. However, using the F-basis of [7], 
a factorized expression for the 2-DA DWPF was obtained. 

In this work, we extend the above result to the TV-DA models, N £ {2,3,4}. 
Our results are restricted to N G {2,3,4} because our proofs require the explicit 
expressions of the weights, while the number of vertices grows ~ 0{N^). However, 
our results are quite simple and have a uniform dependence on N, which allows us 
to conjecture that our expression extends to all N > 2. 

Non- invar iance under state variable conjugation. Our proofs rely on 
the non-invariance of the A^-DA models under state variable conjugation (for non- 
vanishing external fields) . This leads us to look for other models that are similarly 
non-invariant. 
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The sl{r+l\s + l) Perk-Schultz (PS) models, {r, s e N} [8], form another 
class of models that are non-invariant under state variable conjugation, in this case 
because the state variables belong to two different sets with different statistics. 

In these models, the definition of domain wall boundary conditions (DWBC's) is 
not unique. From experience with the A''-DA models, we propose a definition that 
leads to factorized DWPF's. The symmetries of the Perk-Schultz models are such 
that the result is independent of {r,s}. 

Outline of paper. In section 1, we recall basic definitions, introduce the N- 
DA models and obtain the corresponding factorized DWPF. In section 2, we do the 
same for the sl{r + l\s + 1) Perk-Schultz models. Section 3 contains brief remarks 
and appendix A lists the weights of the A''-DA models, for N G {2, 3, 4}. 

1. The A''-state Deguchi-Akutsu {N-DA) models 



1.1. The lattice. We work on a square 
lattice consisting of L vertical and L hor- 
izontal lines, label the vertical lines from 
left to right and the horizontal from top 
to bottom. 

We assign the i-th vertical line an ori- 
entation from bottom to top, a complex 
rapidity variable Ui and a complex ex- 
ternal field variable a^. We assign the 
j-th horizontal line an orientation from 
left to right, a complex rapidity variable 
Vj and a complex external field variable 
A-- 



Vl,(}l 



VL,f3L 



ai 



Ul 



Figure 1. An L x L square lattice, 
with oriented lines and variables. 



1.2. Vertices. Each lattice line inter- 
sects with L other lines. A line segment 
between two intersections is a bond. To 
each bond, we assign a state variable 
a G {1,2,- •• ,N}. The intersection of 
the i-th vertical line and the j-th hori- 
zontal line, together with the four bonds 
adjacent to it, and the set of state vari- 
ables on these bonds, is a vertex Vij. 



v,(3 





t-1 


K2 


t2 







u, a 



Figure 2. The vertex corresponding 
toX^^^{u-v)Xl.- 



1.3. Weights. To each vertex Vij we assign a weight tWjj, that depends on the 

state variables on the four bonds of that vertex, the difference of rapidity variables 
flowing through the vertex, and the two external field variables flowing through 
the vertex. Specifically, a vertex with vertical rapidity and external field variable 
{u, a}, horizontal rapidity and external field variable and state variables 

{ii, (-2, Ki, K2} is assigned the weight Xa,p(^ — v)''^^''^^. These weights satisfy the 
Yang-Baxter equation: 
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J2 X^dn - vr^X^„.,{u - w)ll'XXp,,{v - wf^l'^^l = (1) 

, A2 , A3 

Expressions for all Xa.piu — given in [5], for N e {2,3,4}. For 

completeness we include them in appendix A. From these expressions, one can 
check that the weights of the A^-DA model are not invariant under conjugating the 
state variables. This is due to the presence of external fields {a, /?}. 

Switching off the external fields restores the symmetry of the weights (up to 
global gauge transformations). The (symmetrized) no-field weights coincide with 
the weights of the spin--^^--!- models at their free fermion point. This can be checked 
trivially for A = 2 by setting a = f3 = \/— T. 



1.4. Minimal and maximal state 
variables. We refer to the state vari- 
able cr = 1 as minimal, and to cr = A 
as maximal. 

1.5. The c+ vertex. We refer to the 
unique vertex with minimal state vari- 
ables incoming from the left and ex- 
iting from above, maximal state vari- 
ables incoming from below and exiting 
from the right, as shown in Figure [3l 
as the c+ vertex. In the A-DA models 
c+{a, /3,u-v)= Xa,p{u - w)};^. 





1 


1 


A 




A 



Figure 3. The c+ vertex. 



1.6. Domain wall boundary condi- 
tions (DWBC). We define DWBC's in 
the A-DA model as a form of expanded 
c+ vertex: All (boundary) bonds on the 
left and top carry minimal state vari- 
ables, while all bonds on the right and 
below carry maximal state variables. 



_A 
_A 
_A 
_A 

A 



A A A A A 



Figure 4. N-DA DWBC's. 



1.7. Line-permuting vertices. Any model, on a finite lattice with DWBC's, has 
a pair of vertices that can be used to permute adjacent lattice lines, as we will see 
in detail below. We call these vertices {a+, a_}. In the A-DA models, 

a+{a,l3,u- v) = Xa^piu- v)l'^^ 

(a \ V / \N,N y ) 

a-(a,P,u- V) = Xajj[u ~ w)jv> 
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1.8. Domain wall partition function (DWPF). As always, the DWPF on an 
L X L lattice, Z£^^, is the sum over all weighted configurations that satisfy the 
DWBC. The weight of each configuration is the product of the weights of the 
vertices: 



= E 

configurations 



n 



(3) 



1.9. Properties of the N-BA DWPF. The TV-DA DWPF, N e {2,3,4}, sat- 
isfies four properties: 



Property 1. From the DWBC's, Zg^^ has the form 



ZE^L (W-{/3},W,W) =e(^-i)"^p (W,{/3},W,W) 



(4) 



where p is a polynomial of degree (L — 1)(A'^ — 1) in . 

This can be seen from the fact that the rapidity ui only appears in the left-most 
column, and every vertex in that column is of the form Xa^_/3^.(ui — vjY^'^y From 
appendix A, one can check that 



A 



(5) 



where q[''^ is a polynomial of degree (TV — 1 + <, — A) in . Property 1 then follows 
by noticing that every lattice configuration in the DWPF receives a contribution of 



from the left-most column. 



N 



(6) 



Property 2. Zj^^j^ has zeros in e"i at the {L - 1){N - 1) points 



je{l,...,JV-l},fce{2, 
11111 



N 
N 



X 



_N 
_N 
_N 
_N 

N 



N 



N 



N 



(7) 



Figure 5. Attaching an o_ vertex. 
This can be seen as follows: Multiply .^f^fj;, ^{a}, {/3}, {u}, {u} j with the vertex 

a^{ai,a2, ui — U2)- This corresponds to attaching a type a_ vertex from below, as 
in Figure 5. 

Using the Yang-Baxter equation, we slide the inserted vertex through the lattice 
until it emerges from the top, as in Figure 6. 

The inserted a-{ai,a2, ui — U2) vertex emerges as an o+(ai, 0:2, mi — U2) vertex, 
and in the process, the two left-most vertical lattice lines are permuted. 
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Figure 6. Extracting an a-|_ vertex. 



We conclude that 



zE^L (W,{/3},W,W) 



a+{ai,a2,ui - U2) 



a_(ai, Q!2, Ui — U2) 
ZExL ({"2, ai, • ■ ■}, {/?}, {U2, ui, . . .}, {w}] 

Iterating tlie above procedure {L — 1) times, we obtain 



(8) 



^ \^a^{ai,aj,ui- Uj) 

^L>fL ({"2, ■ • ■ , aL,ai}, {/?}, {U2, . . . , UL,Ui}, {«} j 

The locations of the (L — l)(iV — 1) zeros in e"^ follow from Equation[9] 



(9) 



Property 3. Zj^^j^ obeys the recursion relation 



■'LxL 



N-l N^l 



3 = 1 
L 



n n [^-p'^'PLPkc''--"^] n ( 



3=1 fc=i 



fe=2 



yDA,(lL) 
'(i-l)x(L-l) 

(10) 



where ^'JI'IJ^Jl-i) is the DWPF on an (L— 1) x (L— 1) lattice, with the omission of 
external field variables {ai,/3i} and rapidities {ui,vl\- This is seen by noting the 
lower-left vertex must be Xa^^ij^{ui —vl)i%, which, as can be verified in appendix 
A, satisfies 







0, unless l — I, k — N 



(11) 



Hence, setting e"^ = ^e"^ in ZE^l freezes the lower-left vertex to a type c+, the 
remainder of the bottom row to type a_, and the remainder of the left- most column 
to type a+. Equation 1101 follows from these considerations. 
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Property 4. The DWPF on a 1 x 1 lattice is given by the c+ vertex 



ZDA ^ ^(A^-l)K-..) -Q (12) 
j = l 

which follows from the definition of the DWBC and the weights. 

Lemma 1. The above four properties determine the iV-DA DWPF, N G {2, 3, 4}, 
uniquely. 



Proof. Write Zf^^ for the actual DWPF, and assume there exists some other 
^LxL which satisfies all of the preceding four properties. By Property 4, we have 
Zf'^ = H^^i, which is the basis for induction. Fix an integer n > 2. From 
Properties 1 and 2, n^^„ must be equal to Z^^^^, up to a multiplicative term, C, 
that does not depend on . 

From Property 3 and the inductive assumption ■^(^^i)x(„-i) = ^fnti)x(n-i)' 
we find that the multiplicative constant C = 1. Hence, Z^^^ — n^>f„, proving the 
uniqueness claim by induction. 

1.10. Evaluation of the iV-DA DWPF, N € {2,3,4}. We postulate an ex- 
pression for Z£^j^, then show that it satisfies the four properties of the previous 
section. 



Lemma 2. 





( N-i ^ 




zE^L = n 


^{N-l)j{u,-v,) -Q 
k=l 


-p^-ia2^1-p'=-i/32 

> 


X 


n n 

l<i<j<L k=l 




^] (13) 



Proof. By inspection, the product expression in Equation 1131 satisfies Property 1 
and 4. It contains a factor of Y[f=2 Il^i^ ~ P^~^^'i^3^^^'^^ ] ' which means it 
possesses the (L— ~ 1) zeros required by Property 2. Finally, working directly 
from Equation 1131 we obtain 



L N-\ 

k=\ j=l 

n n (i-/-^/3.fte— ) n (i-/-«ia.e--) ^si^Vi) 

k=l i=l j=2 

(14) 
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N-1 



,j.,^PL,.r = 1 — 1 e ^ - ' "x 



Evaluating Equation [14] at the point e"^ — ^e"^, we obtain 

L N-1 

n (e(^-i)-) n ^w^/T"^ 

fc=2 j=l 
N-lL-1 L 



DA.{1L) 

)x(L-l) 

fe=l 1=1 j=2 

(15) 

Rearranging factors in Equation 1 151 one recovers Equation llOl in Property 3, as 
required. This concludes our proof that the A^-DA DWPF, N £ {2, 3, 4}, factorizes. 

□ 

Conjecture. Equation ll3l is valid for all N > 2. This conjecture is based on the 
fact that our results for N G {2,3,4} have a uniform dependence on N, in which 
only the c+ and line-permuting vertices appear. A study of these very vertices for a 
few values of > 4, indicates that they have analogous forms and properties, hence 
our conjecture. However, a proof of this conjecture requires detailed knowledge of 
all vertices, for all > 4, which is beyond the scope of this work. 

2. The sl{r+l\s + l) Perk-Schultz (PS) models 

This section will be brief, as the arguments are the same as for the A^-DA models. 
The PS models are defined on the same lattice, with the same orientations as the 
A^-DA models. The differences are in the weights, as those of the PS models do 
not depend on external field variables. 

2.1. Two sets of state variables. Following [9], we define two sets i?_={l,. . . ,s+ 
1}, _B+={s + 2,. . . ,r -I- s -I- 2}, and their union 

B = {l,...,s + l,s + 2,...,r + s + 2 = N} 

" V ' " V ' 

2.2. The weights. Let a,b G B. The non- vanishing weights of the sl{r + l|s + 1) 
Perk-Schultz (PS) models are 

jja.af \ J smn?7 ' 

^a,a\^) — \ sinhr?(l+ii) ^ o 

Sinn 77 ' ' 



-Siiili^, a,be B^ OT a,be B+ 

siiih 77n ii 

■ , otherwise 



K>) ^ { ^r;; : < I (16) 

where 77 is a crossing parameter. The labelling of the vertices follows the same 
convention as in Figure 2. 

2.3. Non-invariance under state variable conjugation. It is clear, by inspec- 
tion, that the weights of the sl{r + l\s + 1) PS model are not invariant under 
conjugation of state variables, where A^ = r + s + 2. 
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2.4. Symmetries and {r, s}-independence. From Eauations ll61 it is clear tiiat 
the PS weights arc symmetric in ct- G _B_ and (separately) in the a+ G sjl. One 
can also see that, choosing any (T_ G -B- and any cr+ G as domain wall bound- 
ary variables, no other state variables appear in the domain wall configurations, 
and one obtains a factorized DWPF that is independent of {r, s}, in other words, 
the same result one obtains in the sZ(l|l) modeO. 

2.5. PS Domain wall boundary conditions. We define the sl{r + l|s + 1) PS 
DWBC as follows: The state variables on all bonds on the right and lower bound- 
aries are maximal, a = N = r + s + 2€ i?+ , and the state variables on all bonds 
on the left and upper boundaries are minimal, <t = 1 G -B- . 

Notice that, using the symmetries of the weights, we could have taken any state 
variable a G B- on the left and top, and any state variable a G 5+ on the right 
and below. The advantage of the above choice is that the labels of the c+ and 
line-permuting vertices are precisely the same as those in the A'^-DA models. 

2.6. Properties of the s/(r-|-l|s + l) PS DWPF. The sZ(r -I- l|s + 1) PS DWPF, 
■^LxL' satisfies the following four properties, the proofs of which are precisely anal- 
ogous to those of the A^-DA models. 

Property 1. Given the PS weights and DWBC's, and writing Ui = e''"S Zf^^ 
has the form 



zlL (w, w) = ^r^+'p (m,m) (17) 

where p is a polynomial of degree (i — 1) in [/^. 

Property 2. Using the line-permuting vertices, {R\'\,R^'^}, and the Yang- 
Baxter equations, it is straightforward to show that 



which gives the (L — 1) zeros of p. 

Property 3. Setting ui = vl, we freeze the lower left-hand corner to an R\' 
and obtain the recursion relation 



r^PS I 



'L-1 



' L 



Z 



PS,{1L) 
(L-l)x(L-l) 



Property 4. The initial condition is given by the c+ vertex 

Lemma 3. The above four properties determine the PS DWPF uniquely. 
Proof. The proof is identical to that of Lemma 1. 



(18) 



'^sl{N) Belavin models are analogously symmetric, but all state variables take values in one 
set only ct e {1, ■ • • ,N}. 

^In sl{N) Belavin models, choosing any two distinct state variables to impose DWBC's, no 
other state variables appear in the configurations, and one obtains Izergin's determinant expression 
for the spin-i model, which corresponds to = 2. 
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2.7. Evaluation of the PS DWPF. We postulate an expression for Z^^j^, tlien 
show that it satisfies the four properties of the previous section. 

Lemma 4. 



yPS 



L 



Kk=l 



n ^M("i - Uj)R\\{vj - Vi) 
^ l<i<j<L 



(19) 



Proof. The proof is identical to that of Lemma 2 in the case of the A^-DA models, as 
one can show that the expression in Equation[T9]obeys the required four propertied. 



3. Remarks 

Domain wall partition functions are important in physics because of their role 
in the algebraic Bethe ansatz approach to correlation functions [11], and in com- 
binatorics, because in the special case of the spin-^ model, they lead to counting 
alternating sign matrices [12]. 

The factorizing DWPF's discussed in this work do not lead to new combinatorics: 
One can define corresponding combinatorial objects, but the weights do not lead 
to 1-counting (they cannot all be simultaneously set to 1). The factorization of the 
DWPF's (probably) reflects the fcrmionic nature of the underlying models and the 
absence of interesting counting is in turn a reflection of that nature. 

On the other hand, factorized DWPF's are easier to handle than determinants 
and one expects (from experience with the spin-^ model at the free fermion point) 
that computing the corresponding correlation functions will be easier than in generic 
models. 

The main point of this work is to observe that, if the weights of the line-permuting 
vertices have different zeros, then the corresponding DWPF's factorize. The models 
discussed in this paper offer trigonometric vertex examples of this observation. It is 
possible that all such models are either fcrmionic (with or without external flelds) 
as in the A^-DA models, or contain fcrmions, as in the sl{r + l\s + 1) PS models. 



Appendix A 

In the following, w = u — v, x = e'^ and p = , where {n, N} are co-prime. 
There is no crossing parameter because these models are external field deformations 
of spin- ^ models with a crossing parameter set to the free fermion value. 

N=2. (6 vertices) 



^a,/3(w)};l = (1 - a/Sx) 
Xa^f3{w)l'^^ = (a - Px) 
X^^piw)lj = v/(l-a2) (1-/32) 



Xc.Mw)li = 2:^(1 -«2)(l-/32) 
Xa,f3{w)l'^ = (/3 - ax) 
Xa,f3{w)l'^ ^ {x - afi) 



In [10], a determinant expression for the DWPF of the PS model was obtained. It is 

straightforward to show that it evaluates to the above product expression. 
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N=3. (19 vertices) 

Xa,/3Ml'^ '■ 

Xa.(j{w)l'^ - 

^a,/3(w)i;3 : 

Xa,l3{w)^]2 

N=4. (44 vertices) 

Xa,0{w)\'^ ■■ 

Xa.^0{w)l-ll - 
Xa,(3{'w)lll ■■ 

^a,/3(w)3;2 = 



(1 — a/?x)(l — a(3px) 

0:^(1 -a2)(l-/32)(i _ aPpx) 

{a — /3a;) (1 — appx) 

a;V(l - - - - /3V) 

V(l-a2)(l-/32p) ^/~^' x(a - Px) 

(a — I3x) {a — (3px) 

{j3 — ax){\ — aPpx) 
^(l_a2)(i_/32)(i _ 

V(l-a2p)(l-/32) ^//~^' x(/3 - ax) 

(1 - a^){l - f5^p)x -{[3- ax){(3x - ap) 
^{l-a^){l-P^p) ^lZ^ {a - (3x) 

x{x - a/3)V(l-Q!V)(l-/3V) 
{l + p){a-px){x-ap) 
{(3 — ax){p — apx) 

x/(l-/32)(l-a2p) ^l^ (/3 - ax) 
V(l-a2)(l-a2p)(i_^2)(i_^2^) 

(1 + p)(/3 - a;a)(a; - aP) 
^{l-a^p){l- I3^p){x - a/3) 
(a; — aP){x — aPp) 



(1 - a(3x){l - a(3px){l - app^x) 



a;V(l-Q!2)(l-/32)(l - Q!/3pa;)(l - ajSp^x) 

{a - /3a;) (1 - a/3pa;)(l - a/3p^x) 

a;V(l - - aV)(l - /3^)(1 - /32p)(l - Q;/3p2a;) 



v/(l-a2) (1-/32^) 



x(q! — /3a;) (1 — af3p'^x) 



(a — /3a;) (a — /3pa;)(l — a(3p^x) 

XV(1 - - a»(l - a2p2)y'(l _ ;32)(1 _ ^2p)(l _ ^2^2) 

v/(1 - a2)(l - a2p)(l _ ^2)(l _ ^2^) V^^^2(^ _ 



V(l - a2)(l - /?2p2) 



:x(q — /3x)(a — (}px) 
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Xa,!3{w 



XaA'^ 



Xa,l3{w 
Xa,0{'W 



XaA'W 



Xa,piw 

Xccpiw 

Xa,l3{w 

Xa.piw 
Xa,l}{w 



1,4 
'4,1 
2,1 
1,2 
2,1 
2,1 

,2,2 
1,3 

|2,2 
'2,2 

2,2 
3,1 



,2,3 
1,4 

,2,3 
'2,3 
.2,3 
'3,2 

,2,3 
'4,1 

,2,4 
'2,4 

^,2,4 
^3,3 

,2,4 
14,2 

,3,1 
1,3 

3,1 
'2.2 

,3,1 
3,1 

,3,2 
1,4 

,3,2 
'2,3 
3,2 
'3,2 



(a — (5x){a — Ppx){a — (ip^x) 
{fi — ax){l — a(3px){l — afip^x) 



V(1-Q;2)(1-/?2)(1 - al3px){l - aPp^x) 



- a2p)(l - /j^) ^ / ^ =x{l3 - ax){l - app'^x) 
((1 - a^)(l - P'^p)x -{13- ax){Px - ap))(l - a/Sp'^x) 



V(l-a^)(l-/?» ^/ - (3x){l - aPp^x) 



V(l - - a2p)(l - - P^p)^^^^2={a-Px){l-aPp'^x) 



y/{l - a^p){l - p^p)iil-p^)il-a'^p'^)-{l+p)x{ax~l3p){a-Px)) 
{a - I3x){{l - a^P^){l - p^)x - p{ax - (3p){a - Px)) 



^/ ^V (l-a2)(l-/3V)(a - px){a - ppx) 



= xVl - - Q;V^)(1 - - P'^p'^){x - a/3) 

V(l-p2)(l_p3) 



1-p 



v/(l-aV)(l-/32p2)(a; - a/3) (a - /3a;) 



(1 - p^){x - a/3)(a - Px){a - /3pa;) 
(/? — ax){p — apa;)(l — aPp^x) 



v/(l-/32)(l-a2p) ^/ ^' (/3 - aa;)(l - aPp^x) 

vi-p 

= V(l - Q!^)(l - a2p)(l - /32)(1 - /32p)(l - a/Jp^a;) 

iP - ax){{l - a^P^){l - p^)x - p{P - ax){Px - ap)) 



y(T^^W^^^((l-a') (/3-aa;) (/3a;-ap)) 



XaA'^. 

Xa,(3iw 

Xa^piw 
Xa,0{'W 

XocA'^. 
Xa,f3iw 

Xa^piw 



,3,2 
'4,1 



,3,3 
'2,4 

3,3 
3,3 

3,3 
4,2 

3,4 
3,4 

3,4 
4,3 

4,1 
1,4 

4,1 
2,3 



V(l - n:2)(l - a^p){l - p^p){l - p^p^) ^}-^ {a - Px) 



^ ^(1 P")(^ _g!)^(i_a2p2)(i_^2p)(a. _ c,/3)(a _ f3x) 
1 - p 

{{l-a^p)(l-P^p^)x~{l + p + p^){P-ax)(Px~ap)){x-aP) 
V(l - a2p)(l - p^p) - ap){a - Px) 



x\J 1 — a^p^^/l — P'^p'^{x — a/3) (a; — a/3p) 
1 -P^ 

—{x - ap)ix - aPp)(a - px) 

1- p 

{P — ax){P — apx){P — ap^x) 



V(l - aV)(l - 



p- 



(/3 — ax){a — /3a;) 



12 



^a,/3(w)3;2 
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,3 



v/(l - a^)il - a^p){l - - _ ^x) 

V 1 



1 - p3 



1-p 



-{x — af3)((3 — ax){P — apx) 

V(i-P^)(i-P^), 



J(l-a^p^){l- I3^p)-^ PA _ aP)((3 - ax) 

1 - P 

V(l - «2)(1 - a^p){l - /32p)(l _ /32p2)(^ _ 

1 -p3 



-{x — aP)(x — a(3p){[3 — ax) 



y/{l-a^p^){l~ I3^p^){x - al3){x - a(3p) 
{x — aP){x — a(3p){x ~ a[3p^) 



